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ABSTRACT 

We study the spectral and morphological characteristics of the diffuse Galactic emission in the 
WMAP temperature data using a template-based multi-linear regression, and obtain the following 
results. 1. We confirm previous observations of a bump in the dust-correlated spectrum, consistent 
with the Draine & Lazarian spinning dust model. 2. We also confirm the "haze" signal in the 
inner Galaxy, and argue that it does not follow a free-free spectrum as first thought, but instead is 
synchrotron emission from a hard electron cosmic-ray population. 3. In a departure from previous 
work, we allow the spectrum of Ha-correlated emission (which is used to trace the free-free component) 
to float in the fit, and find that it does not follow the expected free-free spectrum. Instead there is a 
bump near 50 GHz, modifying the spectrum at the 20% level, which we speculate is caused by spinning 
dust in the warm ionized medium. 4. The derived cross-correlation spectra are not sensitive to the 
map zero points, but are sensitive to the choice of CMB estimator. In cases where the CMB estimator 
is derived by minimizing variance of a linear combination of the WMAP bands, we show that a bias 
proportional to the cross-correlation of each template and the true CMB is always present. This bias 
can be larger than any of the foreground signals in some bands. 5. Lastly, we consider the frequency 
coverage and sensitivity of the Planck mission, and suggest linear combination coefficients for the 
CMB template that will reduce both the statistical and systematic uncertainty in the synchrotron 
and haze spectra by more than an order of magnitude. 

Subject headings: diffuse radiation — dust, extinction — ISM: clouds — radiation mechanisms: non- 
thermal — radio continuum: ISM 
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1. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) by the Wilkinson Microwave Anisotropy Probe 
(WMAP) have revolutionized our understanding of cos- 
mology and placed strong const r aints o n cosmological 
param eters (S pe r gel et al.l 120031 120071 : iTegmark et al.l 
2004). Moreover, the WMAP foreground signal rep- 
resents the most detailed and sensitive full-sky maps 
of Galactic microwave emission, providing an enormous 
wealth of information about the physical processes in the 
interstellar medium (ISM). 

1.1. Galactic emission mechanisms 

There are three well established types of Galactic fore- 
ground signals at WMAP frequencies: free-free (or ther- 
mal bremsstrahlung) emission from interaction of free 
electrons with ions, dust emission from grains heated by 
the surrounding radiation field, and synchrotron emis- 
sion from supernova shock accelerated electrons. In ad- 
dition, there are two other emission mechanisms which 
have proven more difficult t o characterize: sp inning; dust 
and the anomalous "haze" (|Finkbeinerl 120041 ) . Spinning 
dust refers to emission from the smallest dust grains 
which have non-negligible electric dipole moments and 
are excited into rotational modes through a variety of 
mechanisms (cf. iDraine fc Lazarianl[l9 98b). The physi- 
cal origin of the haze is uncertain. 

Though initially controversial ([Bennett et al.l 120031 ). 
numerous authors have presented evidence for a spin- 
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ning dust spectrum wh en combining WMAP data 
with external data sets (|de Oliyeira-Costa et al.l [ 2004: 
Finkbeiner. Langston. fc Minted [2004: Bough njfc Poberl 



200 ?]). Using only WMA P data, iBennett et all (|2i / 
and iHinshaw et all {2007) point out that it is difficult 
to spectrally distinguish certain spinning dust models 
from synchrotron. Howeve r, in our companion paper 
(jDobler fc Finkbeinerl [20071 . hereafter DF07), we show 
that a spinning dust spectrum is indeed recoverable us- 
ing exclusively WMAP data, though it is not spatially 

correlated with the thermal dust e mission. 

The haze was originally thought (|Finkbeinerl l2004h to 
be free-free emission from ionized gas which is too hot 
to be visible in recombination line maps and too cold to 
be visible in X-ray maps. However, gas at the r equired 
temp erature T ~ 10 5 K is thermally unstable (|Spitzerl 
Il978f ). Furthermore, we will show in §3]that the spectrum 
of the haze is inconsistent with free-free emission and is 
most likely explained as a hard synchrotron component 
which is morphologically and spectrally distinct from the 
above mentioned soft synchrotron. 

1.2. MEM analysis 

With the first and third year data releases, the WMAP 
team provided a dual foreground analysis: a maximum 
entropy method (MEM) and a template fitting algo- 
rithm. The former was intended to improve our un- 
derstanding of the astrophysics of foreground emission 
while the latter (more statistically stable) algorithm was 
meant to produce CMB maps with well characterized 
noise properties for use in the cosmolo gical analysis 
([Bennett et al.ll2003l : IHinshaw et al.ll2QQ7[ ). 

Although the total observed emission matches the 
MEM model to better than 1%, "Low residual solu- 
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tions are hig hly constrained, but not necessarily unique 
or correct" ([Bennett et al.l I2QQ3L p. 108). The MEM 
method is a pixel-by-pixel fit w hich minimizes th e func- 
tional H(p) = A(p) + \{p)B{p) (jPress et al.lll992h where 
A is a regularizing parameter, 

A(p) = X \p) = YF(y,P) ~ T m (v,p)] 2 /a 2 , (1) 

V 

B(p) = J2 T c (p) ln[T c (p)/P c (p)}, (2) 

C 

and the sum is over Galactic emission components. 
P c (p) is a prior template for com ponent c, normalized 
to the same frequen cy as T c (see iBennett et al.l l2QQ3t 
iHinshaw et alJl2QQ7l for details). 

The noise properties of a MEM-derived CMB map are 
complicated, e.g. by the fact that noise is clipped to 
be non-negative by the logarithm in Equation [2l Be- 
cause simple noise properties are desirable for a cosmo- 
logical power spectrum analy sis, a template - based map 
is used instead . In t he end, IBennett et al.l ([2003) and 
IHinshaw et al.l ([2007) fit spatial templates for the dust, 
free-free, and synchrotron ISM emissions to WMAP Q, 
V, and W bands. They find that the remaining contam- 
ination in the maps is sufficiently subdominant outside 
the Kp2 mask so that the effect on the CMB power spec- 
trum is negligible. 

1.3. "Anomalous" ISM emission mechanisms 

Given their dual approach, an important question to 
address is: why have the foreground analyses of the 
WMAP team not positively identified the sp inning dust 
emissi on in DF07 and the haze emission in Finkbeine3 
([2001 ? 

With the MEM analysis, there are two pitfalls. First 
and most importantly, the MEM analysis generates a 
map of the spectral behavior of the soft synchrotron in 
each pixel Since the free-free and dust spectra are kept 
fixed in the model, deviations away from these assumed 
prior spectra are absorbed into this "spectral index map" 
for synchrotron. Second, when minimizing H(p), if the 
recovered model T m (z/,p) yields negative pixel values, 
X(p) is increased until the results are greater than zero. 
These two effects combine so that foreground emission 
which does not match the prior templates and spectra 
(i.e., spinning dust and haze emission) are simultane- 
ously absorbed into the synchrotron spectral index map 
and washed out by adding priors back into each pixel to 
enforce positivity. Since minimizing H(p) in this way is 
repeated multiple times (with the synchrotron spectral 
index map being "updated" with each iteration) the re- 
sults of the MEM analysis naturally strongly resemble 
the priors. 

With the template fitting algorithm, the difficulty 
arises in the choice of a synchrotron template. Because 
the WMAP sign al to noise ra t io is far superior to pre- 
vious surveys . IHinshaw et al.l (|2007l ) use the difference 
of the two lowest WMAP frequency maps (K— Ka) as a 
template for synchrotron 3 . They acknowledge that this 
template necessarily contains free-free emission as well, 

3 IBennett et al. ( 2003) used the synchrotron template described 
in £12.11 and an excess south of the Galactic Center is indeed visible 
in their residual maps (their Fig. 11, upper right panel). 



but it also contains the haze and spinning dust emission. 
Thus, the haze and spinning dust are simultaneously ex- 
plicitly fit with this template. In addition, the spectra of 
the dust emission and free-free emission are again kept 
fixed so deviations are difficult to identify. 

Given the numerous challenges involved in a CMB fore- 
ground analysis and the serious questions about whether 
all the relevant foreground emission mechanisms have 
even been identified yet, we choose an approach that is 
simple enough to have well characterized noise proper- 
ties, but flexible enough to allow us to find surprises. 
In the limit where the spectrum of each component is 
invariant with position, one still has the choice of assum- 
ing a perfect spatial template, or assuming knowledge of 
the spectrum of each component. Given that the spec- 
trum of each component can vary with position, neither 
of these approaches is strictly correct. Nevertheless, as 
we shall see, there is still much to be learned by making 
a too rigid assumption and then studying the resulting 
residuals. 

2. METHODS 

We proceed under the assumption that the morpholo- 
gies of the foreground components are well characterized 
by external data sets and that their spectra do not vary 
significantly over our regions of interest (see §3.2p . Thus, 
in principle, we should be able to find the appropriate lin- 
ear combination of the foreground templates which, when 
combined with a CMB template, can be subtracted from 
the WMAP data leaving only random residuals consis- 
tent with noise. 

2.1. Template maps 

Each of the templates used in o ur fits has been dis- 
cussed in detail by previous authors ([Be nnett e t al .112003 ; 
lFinkbeinerll2004l : (Davies et aD2QQd IHinshaw et alJ2QQ7h 
and so we only briefly review them here. 

Free-free: Free-free (or thermal bremsstrahlung) emis- 
sion originates from the Coulomb interaction of free elec- 
trons with ions in a warm gas. Since this emission is 
proportional to gas density squared, maps of Ha re- 
combination line emission (which is also proportional to 
density squared) roughly trace the morphology of the 
gas and thus also the free-free emission. Our template 
for this foreground is the Ha map, assem bled from the 
VTS S (jDennison et"aIj[T998h. SH ASSA (paustad et all 
2001 ), and WHA M (jHaffner et all HIl) surveys by 
Finkbeinerl (|2003l ). Extinction by dust (both in front 
of and mixed with the warm gas) presents a potential 
chall enge in interpreting Ha as a tracer of fr ee-free emis- 
sion (Bennet t et al.ll2003l : iFinkbeinerl 12004). We correct 
th e Ha map for dust extinction using the prescription 
in IFinkbeinerl (|2003l ). and further mitigate the effect by 
limiting our interest only to regions where the dust ex- 
tinction is A(Ha) = 2.65E(B — V) < 1 mag, where A is 
related to the dust optical depth = A/(1.086 mag). 

The free-free spectrum is co nstrained by physics 
(|Spitzerlll978l : IBennett et aLll20Q3h as 

Tocz/*, (3) 

where a ~ —2.15, T is in antenna temperature, and the 
proportionality depends only on the electron tempera- 
ture T e . Thus the only free parameter is a scaling factor 
equivalent to the electron temperature on the sky (but 
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see ^3]) , though this temperature is expected to vary with 
position. 

Thermal and spinning dust: The emission produced 
from tiny interstellar dust grains vibrating in equilibrium 
with the surrounding radiation f i eld ha s been mapped 
across the sky by ISchlegel et al.l (ll99cf). We use their 
full- sk y map evaluated at 94 GHz by [Finkbeiner et al.l 
(1999|) (hereafter, the FDS map) as a dust template. 4 
The smallest of these dust grains are expected to have 
a non-negligible electric dipole moment and so can also 
emit radiation at WMAP frequencies through rotational 
modes excited by collisions with ions. Thus our template 
also traces spinning dust emission. Since the spectral 
dependence of the spinning dust is not well known, we 
use our fit to constrain the frequency dependence of the 
dust-correlated foregrounds. 

Soft synchrotron: As relativistic, shock accelerated 
electrons travel through the Galactic magnetic field, 
they emit synchrotron radiation with a characteristic 
frequency dependence oc v& (in antenna temperature 
units). At 4 8 MH z, this emission was measured by 
lHaslam et al.l ([1982), and we use their full-sky map 
as a tracer o f soft synchrot ron. As pointe d out by 
iBennett et "all (|2003f ) and lHinshaw et all (|2007f ) the spec- 
tral index (3 is expected to vary across the sky, and in 
particular, the spectrum may be harder near regions of 
recent supernova activity. Though we use our fit to eval- 
uate the spectra of the 408 MHz-correlated synchrotron 
emission, we note that a value of (3 = —3.05 removes 
most of the emission at high latitude (most notably the 
prominent "North Galactic Spur" feature). 

CMB: As we shall see in ^3 the choice of CMB es- 
timator can dramatically affect the inferred foreground 
spectra in a given fit. To illustrate this sensitivity, we 
use six different CMB estimators defined as follows. 

• CMB1 - The published internal linear combination 
(ILC) map derived by the WMAP team for the 3-yr 
data. 5 

• CMB2 - An ILC with the coefficients that the 
WMAP team have found best canc el their Region 
fore grounds in the three year data lHinshaw et ail 
(|2QQ7T ). This ILC is given by, 

CMB2 = 0.156 Tk - 0.888 T Ka + 0.030 Tq ( . 

+ 2.045 T v - 0.342 T w , W 
where Tj is the observed WMAP temperature data 
in band j, in thermodynamic AT units. 

• CMB3 - An ILC which minimizes the variance over 
our unmasked pixels. Here 

CMB3 = -0.032 T K - 0.205 T Ka + 0.037 Tq 
+ 0.441 T v + 0.760 T w . 



(5) 



• CMB4 - A "high frequency (HF) estimator" that 
removes the dominant foregrounds (thermal dust 
and free-free) from the 94 GHz WMAP data, 

CMB4 = T w — FDS — AHa, (6) 
where the constant A is determined from the 
approximate free-free amplitude at 23 GHz and 
traced to 94 GHz via Equation [3l 

4 the Ha and FD S maps can be found online at 
|http:/ /www , skymaps. info 

available at http://lambda.gsfc.nasa.gov/ 



• CMB5 - A model for the thermal dust is sub- 
tracted from all of the WMAP bands using Td us t = 
(z//94 GHz) 17 x FDS. A minimum variance ILC 
is then generated from this thermal dust pre- 
subtracted data (denoted by primes). The variance 
is minimized for, 

CMB5 



0.104 T^ - 
+ 0.317 



0.289 T^ a - 
+ 1.059 T^, 



0.190 T^ (7) 



• CMB6 - A cleaned map of the 3 -year data, cleaned 
with the iTegmark et "all (|2QQ3f ) (TOH) method. 6 
This method utilizes a linear weighting of the data 
in which the weights depend on the multipole i of 
the spherical harmonic expansion of each of the five 
WMAP bands. 

In the limit where the noise in each WMAP band is 
equal to <jq, the measurement noise of the ILC is simply 



(8) 



so CMB3 and CMB5 are significantly less noisy than 
CMB1 and CMB2. Due to the complicated weighting of 
the TOH method, the measurement noise properties of 
CMB6 are quite complicated. 

Mask: In addition to masking out all point sources 
listed in the WMAP team's three year catalog, as noted 
above, we mask all regions of the sky where the Ha ex- 
tinction due to dust A(Ha) = 2.65E(B — V) > 1 mag. 
We also mask out the LMC, SMC, M31, Orion-Barnard's 
Loop, NGC 5090, and C"Oph. This mask covers 21.5% 
of the sky. 

2.2. Fitting procedure 

Our model is that the observed WMAP data is a linear 
combination of the foreground templates plus the CMB 7 
plus noise. Therefore, we want to solve the matrix equa- 
tion 

Pa = w, (9) 

where w is the CMB-subtracted WMAP data and P is 
a "template matrix" whose columns consist of the fore- 
ground templates outlined in §2.1[ for the coefficient vec- 
tor a whose entries represent the weights of the individual 
foregrounds. 

The template matrix can be represented schematically 
by a block-diagonal matrix, 

//': 

Pi 

P 3 

Pi 



P = 



\ 



(10) 



V 



p 5 J 



where the 5 blocks correspond to the 5 WMAP bands, 
and each block has the form 

( fl,b dl,b $l,b \ 
f2,b d 2 ,b 5 2 ,6 
f3,b d<3,b $3,6 

(11) 



V fN p ,b dN p ,b SN p ,b ) 

6 available at http://space.mit.edu/home/tegmark/wmap.html 

7 Variations in the foreground spectra from place to place on the 
sky will be explored by fitting smaller regions in upcoming sections 
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where f, d, and s are the templates for free-free, ther- 
mal and spinning dust, and soft synchrotron emission 
(in thermodynamic mK) respectively. This makes P a 
bN p x 15 matrix. For each template, the mean of the 
unmasked pixels is subtracted, making the results of this 
fit insensitive to zero-point errors in the templates. For 
each template column, the first index represents the pixel 
number and the second represents a WMAP frequency 
band — i.e., 1=23 GHz (K), 2=33 GHz (Ka), etc. The 
total number of un-masked pixels in each map is N p . For 
our most general fits, we assume no knowledge of the f, 
d, and s spectra, and so those templates do not differ 
for each band. Rather, we explicitly fit the spectra as 
discussed below. 

The CMB-subtracted WMAP data and the coefficient 
vector are column vectors, 



TABLE 1 



/ 



^2,1 

Ts,i 



7i, 



Cl 

C3 



Cl 



\ 



and 



ad,i 

0>d,2 

a s ,2 



\ a «,5 / 



(12) 



where Tij is the observed WMAP temperature data in 
pixel i and band j, and c is one of our CMB estimators. 
Their lengths are 5N p and 15 respectively. 

Since P is not a square matrix with linearly indepen- 
dent rows, it is not invertible. To solve Equation[9]we cal- 
culate P + , where + denotes the pseudoinverse. 8 The so- 
lution a = P + w minimizes the quantity e 2 = | |Pa— w| | 2 , 
so that if we divide both sides of Equation [9] by the un- 
certainty a (following Bennett et al. 2003 and Hinshaw 
et al. 2007, we use the mean measurement noise in each 
WMAP band), the solution 9 



a=(P/a) + (w/a) 
minimizes the quantity 

HPa - wll 



p 

— a — 


w 


a 


a 



(13) 



(14) 



This fitting procedure is flexible in that additional 
foreground components can be incorporated by simply 
adding columns to the template matrix. We will exploit 
this feature in ^3] to model the anomalous "haze" excess 
emission towards the Galactic center. 

The five components of each coefficient (e.g., a^j, with 
j= [1:5]) represent a fit of both the amplitude and spec- 
trum of the associated emission. By simultaneously fit- 
ting all of the spectra for all of the foregrounds, we can 

8 The pseudoinverse is denned as, P + = £>X+?7 T , where the sin- 
gular value decomposition of P = UJ2B T and E+ is the transpose 
of X with all non-zero singular values replaced by their inverse. In 
the case of a square, non-singular matrix, P + — ► P _1 . 

9 From the properties of the pseudoinverse, PP+P = P and 
[PP+] T = PP+, it is easy to show that, if the columns of P are lin- 
early independent as they are in our case, P + = [P T P] _1 P T , mak- 
ing our technique equivalent to other x 2 minimization t echniques 
(e.g. ITegmarke t al. 2003: Ide Oliveira-Costa fe Teg mark 20(H). 
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Note. — The dif feren t types of fits performed with the 
procedure outlined in ^2.21 using CMB5 for a CMB estimator. 
Fit types are characterized by which spectra are fit, and in 
which region of the sky: full-sky (FS), Galactic center (GC, 
l,b = [—45 : +45]), or specific regions of intere st (R G, see 
Figure [3]). f The haze template is described in £13.11 J The 
X 2 /v for the RG fits are averaged over all 12 regions, 
completely decouple the bands from each other in our 
fits. Additionally, we also perform less general fits in 
which various spectra are fixed to follow the dependen- 
cies in §2.11 (see Table [TJ. We point out that, fitting the 
spectrum of the free-free emission serves as a check on 
the assumption that it is well described by Equation [3] 
(see g3j). 

3. RESULTS 

We have performed both full-sky fits as well as fits of 
smaller regions — the motivation being that both the soft 
synchrotron and spinning dust spectra should vary from 
place to place across the sky. Our fits are characterized 
by the x 2 statistic of Equation [TH and by the residual 
map, 

r = Pa-w, (15) 

which guides intuition and serves as a visual aid in eval- 
uating the goodness of fit. 

The fits can be separated into eight types as outlined 
in Table [H The types designate which spectra are fit and 
whether the fit was over the full-sky (FS), the Galactic 
center (GC - Z, b = [-45 : +45]), or the "regions" (RG), 
which were chosen as regions of interest based on r for a 
full-sky fit (see Figure [3]). Also shown in Table [1] are the 
X 2 jv for each fit type, where v is the number of degrees 
of freedom and is given by 



N p - N a , 



(16) 



with 7V a the number of elements in the coefficient vector. 
Throughout the rest of this paper, we concentrate exclu- 
sively on our most general fit types 7 and 8, and unless 
otherwise noted, use CMB5 for a CMB estimator. This 
estimator is attractive both because of its uncomplicated 
noise properties and because the thermal dust emission 
has arguably the most well understood amplitude and 
spectral behavior making it amenable to pre-subtraction 
before forming an ILC. 

3.1. Full-sky and Galactic center residual maps 

Figure [1] shows rps7 as well as the unsubtracted maps 
for K, Ka, and Q bands, stretched to ±0.25, ±0.12, and 
±0.08 mK respectively. The fit yields Xfs7/ v = 
removing 95.8%, 95.7%, 96.3%, 97.5%, and 99.7% (for K, 
Ka, Q, V, and W 10 bands respectively) of the variance 

10 The noise in W band is especially low because the coefficient 
of T\y in the CMB5 estimator is very close to one, so the CMB- 
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Fig. 1. — Residual maps r for three of the five bands and for FS7 and FS8 fits using CMB5 (the K, Ka, and Q band maps are 
stretched to ±0.25, ±0.12, and ±0.08 mK respectively; all maps are mean subtracted). The unsubtracted WMAP data are also shown 
for comparison. The FS7 fit removes much of the emission, however there is a remaining excess towards the Galactic center. This excess 
emission is particularly notable south of the Galactic center where obscuration by dust and gas is negligible. FS8 incorporates a simple 
spatial template for this haze and removes much of the residuals in that region. 

from the WMAP data. However, it is clear from the 
second row of Figure [T] that there is still a remaining 
emission residual towards the Galactic center (GC). This 
res idual is the "haze" present in the 1 year data as shown 
bv lFinkbeinerl (|2004ft . 

Although the average power in the haze is small (with 
a mean of just 0.59 kJy/sr per pixel at 23 GHz within 30 
degrees of the GC in the southern sky), our fit may be 
compensating for its presence by adjusting the weights 
of the other foreground templates. To relax the stress on 
the fit, we adopt a crude model for the haze emission, 



h oc 







j_ 

^0 



for r < r ; 
for r > 7*0, 



(17) 



where r is the distance to the Galactic center and we 
arbitrarily set ro = 45 degrees. Since the emission mech- 
anism is unknown, we fit the spectrum of the haze as well. 
The extent to which the other fit parameters change gives 
an idea of the cross talk between the haze and the other 
templates. 

The third row of Figure [T] shows residual maps for FS8. 
It is clear from rpss alone that the fit is improved, partic- 
ularly in the southern GC where obscuration from dust 
and gas is minimal and in the high latitude north where 
the North Galactic Spur synchrotron feature is no longer 
over subtracted at 23 GHz (compare rows 2 and 3 of 
Figured]). Though xV^fss = 2.97 is only slightly lower 
than the FS7 fit, we emphasize that the number of de- 
grees of freedom is quite large (y ~ 155, 000) and so the 
likelihood for the FS8 model is significantly higher. Fur- 
thermore, we are including many pixels at large Galactic 
latitudes where the signal to noise is very low and the 
amplitude of the haze is very small. 

Since the ISM may have different properties near the 
GC (i.e., lower gas temperature due to more efficient 

subtracted W band data has essentially no W band data in it. 
Therefore results from W band should be viewed with suspicion. 




Fig. 2. — The same as Figure [T] except for the GC fits (I, b = 
[—45 : 45] degrees). Though GC7 explicitly fits the free-free, dust, 
and synchrotron spectra, the haze is still present in rQC7 indicating 
that it is morphologically dissimilar to the those templates. The 
quality of the fit is substantially improved with the inclusion of a 
haze template (GC8). The extended structure in the north-west 
GC is likely due to an imperfect haze template. 



cooling because of higher metallicity, increased super- 
nova activity, etc.), it is instructive to consider this re- 
gion separately from the rest of the sky. Figure [2] shows 
r maps for the GC7 and GC8 fits. There is a significant 
decrease in x 2 l v from GC7 to GC8 with the inclusion 
of our haze template — from 6.61 for GC7 to 5.08 for 
GC8. The substantially oversubtracted regions near the 
edge of the mask in rcc7 indicate that the fit is indeed 
attempting to compensate for the haze by adjusting the 
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Fig. 3. — Illustration of the 12 regions independently fit in the 
RG fits. Breaking the sky up in this way removes many of the over- 
and under-subtracted features in the residual maps, particularly 
the excess emission in the north-west GC. 

amplitudes of other templates. Although there is still 
some over subtraction in tqcs , the overall quality of the 
fit is much improved, particularly in the southern sky. 

In the northern sky, there is a large structure just 
north-west of the GC in and around the region of Rho- 
Oph. Typically, extended regions of over- or under- 
subtraction in our fits are indicative of variations in the 
physical conditions of the emission media. However, in 
this circumstance, it is unclear if the under-subtraction 
is due to an over simplification of the radial haze profile 
given our template for this component (Equation [TTj) . 
The next step is to subdivide the sky into certain re- 
gions of interest where our fit is the least successful and 
fit those regions explicitly. 

3.2. Regional fits and composite maps 

Our regions of interest were identified as regions of 
particularly notable over- or under-subtraction in rpss- 
Figure [3] shows the boundaries of these regions superim- 
posed on the raw WMAP data at 23 GHz, rpss, as well 
as a composite map for trgs (a fit which includes our 
haze template). The residual map (which is actually the 
residual maps of each region stitched together with no 
smoothing) is shown for both RG8 and RG7 fits in Fig- 
ure [H (the maps stretched to the same units as Figure 

At all frequencies, the RG8 fit more effectively removes 
the foregrounds than the RG7 fit. There are more re- 



gions of over-subtraction in trg7 (both at high latitudes 
and around the mask edges) compared to trgs, the re- 
gion boundaries are somewhat less continuous in rRG7, 
and most importantly, the haze is still present in the 
trg7 maps. These composite r maps are completely un- 
smoothed and it is a testament to the quality of our RG8 
fit that there are no discernible large scale brightness gra- 
dients between adjacent regions. 
Lastly, we define the residual haze map 

rjy = r RG 8 + a^h, (18) 

where the appropriate amount of haze is added back in to 
each region (i.e., the same amount as was subtracted by 
the fit). The bottom row of Figure [H shows r#. Despite 
the fact that the haze fit coefficients are not constrained 
to be continuous across the region boundaries, r# has no 
clear discontinuities. 

3.3. Foreground spectra 

Figure [5] shows the FS8 spectra (coefficient vectors) 
in kJy/sr per template unit: Rayleighs for Ha, mK for 
FDS and Haslam, and arbitrary units for the haze. The 
most striking feature of these spectra are that they are 
very sensitively dependent on the estimator used for the 
CMB. For example, the synchrotron spectrum actually 
appears to turn up at high frequencies for CMB1, CMB2, 
and CMB3, while CMB4 and CMB5 give the more phys- 
ically motivated power law type spectra. This depen- 
dence is entirely due to the contamination of the CMB 
by foregrounds, which can never fully be removed for 
any CMB estimator. Thus, when we remove the esti- 
mator from the WMAP data to perform the foreground 
fit, we have inadvertently added (or subtracted) some 
foregrounds from the data with essentially the spectrum 



of the CMB (I v 



J ). We emphasize that the con- 



tamination is small, so that it has minimal effect on the 
variance of the CMB estimator, but it is large compared 
to the relative amplitudes of the foregrounds. This is es- 
pecially true for estimators which minimize the variance 
of an ILC; in this case the contamination is proportional 
to the cross correlation of the true CMB with the true 
foregrounds (see Q). 

Despite the large uncertainties, there are concrete con- 
clusions that can be drawn from Figure O First, the 
Ha-correlated emission does not follow the I v oc 
frequency dependence as expected. Instead there is a 
bump in t he spectrum around 30 GH z. In our compan- 
ion paper (jDobler fc Finkbei ner 2007) we argue that the 
Ha-correlated emission has a spectrum that is consistent 



with a classical v 



-0.15 



spectrum plus a WIM spinning 



dust component. Second, although the soft synchrotron 
and haze spectra vary substantially with CMB estimator 
type, for a given type, the haze is always harder than 
the normal soft synchrotron, a point which we explore in 

more detail be low. 

Finally, lik e I Bennett et all (120031), B mkbeinerl (|2004f ). 
iDavies et all (|2QQ6h . and iHinshaw et all ([2007), we find 
that the dust-correlated emission falls from 94 to 61 GHz 
but then rises to 23 GHz consistent with emission from 
both thermal and spinning dust. Since these spectra are 
the result of fits over large areas of the sky and the spin- 
ning dust spectrum is expected to vary with position, it 
is not surprising that we do not see a peak in the dust- 
correlated emission in the range 20-40 GHz as in the 
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Fig. 4. — Residual maps for the RG7 and RG8 fits using CMB5, with the same stretch as Figure ^ The haze is still present in the r*RG7 
fits despite the relatively small size of the fitting regions. The bottom row shows trgs with the subtracted haze in each region added back 
in to that region. These maps are unsmoothed and no continuity constraints are placed on the region boundaries. 



iDraine fc Lazarianl (|l998bl ) models. Rather, we are see- 
ing a superposition of many spinning dust spectra with 
varying peak frequencies. This has led to the misidentifi- 
cation of the this dust-correlated emission as synchrotr on 
in the past ([Bennett et al.ll2QQ3l : iHinshaw et alJl2QQ7h , 

It is tempting to conclude that the RG fits for the indi- 
vidual regions can be used to construct a map of the vari- 
ation in the spinning dust spectra across the sky. Such 
a map would represent an "excitation map" or "irradi- 
ation map" for the dust grains. However, this is simply 
not possible given the contamination of the CMB esti- 
mator by the foregrounds. 11 



3.4. Comparison of the haze and soft synchrotron 

From Figure [5] it seems that the haze spectrum is in- 
consistent with a ^~ 015 free-free spectrum. While some 
of the CMB estimators do yield such hard spectra at low 
frequencies, those fits turn up at high frequencies. This 
behavior is much less likely than simply a contamination 
of the CMB estimator by the haze emission. Further- 
more, the haze is most clearly visible in the southern 
GC where dust obscuration is minimal so that the haze 
should either show up in the Ra map (if it is T < 10 5 
K gas) or in the ROSAT X-ray data (if it is T > 10 6 
K gas). Since neither of these is the case, and since gas 
at intermedia te temperatures is thermally unstable (see 
ISpitzerl fl978). we concluded that the most likely source 
of the haze is synchrotron emission from relativistic elec- 
trons. Figure [5] also shows that this synchrotron emission 
is harder than the soft synchrotron traced by the Haslam 
408 MHz map. 12 

11 Uncertain zero point variation across the maps (due to in- 
strumental limitations, imperfect zodiacal light subtraction, etc.) 
is also a consideration, but it is subdominant. 

12 In their first and third year temperature analysis, the WMAP 
team concludes that the spectrum of the soft synchrotron is harder 
near the Galactic plane and p articularly near the Galactic center 
(Bennett et al. 2003; Hinshaw et al. 2007). This may be inconsis- 
tent with [Page et al. (|2007h and IKogut et all (gooj) who find a 



TABLE 2 



CMB 




Ps 






Ph 




estimator 


23/33 


33/41 


23/41 


23/33 


33/41 


23/41 


1 


-2.86 


-2.58 


-2.81 


-2.14 


-2.13 


-2.23 


2 


-2.78 


-2.47 


-2.71 


-1.98 


-1.98 


-2.01 


3 


-2.75 


-2.53 


-2.66 


-2.24 


-2.56 


-2.31 


4 


-3.15 


-3.33 


-3.22 


-2.57 


-3.14 


-2.76 


5 


-3.10 


-3.13 


-3.15 


-2.39 


-2.67 


-2.52 


6 


-3.00 


-3.05 


-3.01 


-2.37 


-2.92 


-2.51 



Note. — The best fit spectral slopes for th+s (Ps) an d 
r H (Ph) maps for the RG8 fits at 23, 33, and 41 GHz over 
unmasked pixels (with I = [—25 : +25] and b = [—45 : 0], see 
Figure [6j. Though there is significant scatter in the inferred 
spectral index, the haze residual is always harder than the total 
synchrotron signal. 

Figure [6] shows density contours for a scatter plot of the 
unmasked pixel values in r# and the total synchrotron 
residual map, th+s — r H + a s s, in antenna temperature 
for various frequency combinations. The pixels shown 
have / = [-25 : +25] and b = [-45 : 0]. Despite the large 
scatter, the r# emission appears to be a distinctly sepa- 
rate component of synchrotron emission with a spectral 
index that is significantly harder than the synchrotron 
(as shown in Table [2] this behavior persists for all CMB 
estimator types, though the precise spectral indices are 
very uncertain due to the CMB estimator bias describe 
in Q). This point is underscored in an RGB color coded 
map of the 23, 33, and 41 GHz r# and maps. The 
haze emission is distinctly bluer (harder spectral index) 
than the total, redder (softer) synchrotron emission in 
the Galactic center. 

Figure [8] shows the total intensity of the haze as a func- 
tion of distance south of the Galactic center. The radial 
bins are 20 degrees wide and separated by 1 degree in 
latitude. Our simple 1/r profile for the haze emission is 

more or less constant synchrotron spectral index across the sky 
from the WMAP polarization data. Giv en th e difficulties in the 
interpretation of the MEM analysis (see fll.3|) , the reason for the 
discrepancy is unclear. 
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Fig. 5. Foreground spectra for the FS8 fits with our six CMB estimators: solid = CMB1, dotted = CMB2, short dashed = CMB3, 
dot-dashed = CMB4, dot-dot-dot-dashed = CMB5, and long dashed = CMB6. The error bars on the coefficients are the formal error bars 
on the fit (see text). 



not adequately describing the structure of the haze — 
which is also clear from the residual maps in Figures [1] 
and [2J Furthermore, the systematic error bars due to 
chance correlation between the haze and the CMB are 
quite large, particularly at high frequencies. 

4. CMB ESTIMATOR BIAS 

While the statistical uncertainty in our fits comes from 
measurement noise in the WMAP data, systematic un- 
certainties in our fits are dominated by contamination of 
the CMB estimator by the foreground components. We 
have used six different CMB estimators, but they can 
be separated into two categories: ILC type estimators in 
which multiple WMAP bands are combined to approx- 
imately cancel the foregrounds (CMB1, CMB2, CMB3, 
and CMB5), and HF estimators in which a model for the 
thermal dust and free-free emission is removed from the 
highest WMAP band to approximately leave the CMB 



(CMB4). 



4.1. ILC estimators 

iHinshaw etafl (|2007f ) point out that the WMAP data 
consists of the "true" CMB, T c plus some additional 
contamination by foreground components Tj which have 
some spectral dependence Nij where i and j represent 
the observing band and foreground component respec- 
tively. That is, 



Wi 



(19) 



9 




Fig. 6. — Number density contours of unmasked pixel (with / = [—25 : +25] and b = [—45 : 0]) temperatures in the yh (blue) and th+s 
(dashed) maps for the RG8 fits at 23, 33, and 41 GHz for CMB5. Though the spectral slope is somewhat uncertain between each band 
and can vary significantly depending on which CMB estimator is used (see Table [2}, the best fit spectral slope for the haze emission is 
distinctly harder than for the total synchrotron. 
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Fig. 7. — An RGB representation of yh and yh+s for RG8 with CMB5. The color coding indicates the spectral index, in antenna 
temperature, of a given pixel. In particular, the bluer haze region (box) indicates a harder spectrum than the redder synchrotron emission. 
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so that an internal linear combination of the WMAP 
bands is, 



(20) 



where Q are constrained to sum to one in order to pre- 
serve unity response to the CMB, and Tj parameterizes 
the contamination of L by foreground j. 

Let us consider that there are four types of foreground 
emission: free-free (T), dust (D), soft synchrotron (S), 
and the haze (H). Each of these has an associated T, 
and so the ILC is given by 

L = T c + Y F T F + T D T D + T S T S + V H T H . (21) 

Since Q are chosen to minimize the variance in T, 

d{L 2 ) 



= o 



(22) 



for all j, and we have explicitly assumed mean subtracted 
maps so that (L) (T c ) = (Tj) 0. Now, 

T 2 = i? + 2E,-r j r c r j 



+ 2 Ej Efc^j FjTkTkTj 

+ v . r 2 r 2 

' ^ J 3 3 



(23) 



(24) 



so that Equation [221 reads 

= 2<T c T i > + 2^^^) + 2r j (T 2 ). 

(note: the notation in Equation [Ml can be compactified 
into = E^ r ^( T i T ^) where ^ = [c,F,D,S,H] and 
r c = 1) Thus, we can solve explicitly for 



(T c T J ) + Y. k ^k{T J T t 



(Tj) 



(25) 



In the limit of only one foreground Tf the second term 
(which represents cross correlation between the different 
foreground morphologies) disappears and this reduces 
to T = r(T c ,T/) = —a c f/(jf as derived with a similar 
method by Hinshaw et al (2006). 

According to Equation [25J the ILC map is biased 
towards anti-correlation between different "true" fore- 
ground emission morphologies and the "true" CMB. 
However, since we do not know the true T c or Tj a pri- 
ori, we have no information about how much of each 
foreground to add back into L in order to correct for this 
factor. 

We can use Monte Carlo techniques to estimate the 
amplitude (but not the sign) of the error in L due to 
foreground contamination as follows. We construct 100 
realizations of the CMB by generating random phases 
for the measured power, Q, in each Fourier mode of the 
binn ed three year WMAP power spectrum (^ max = 986) 
from ISpergel et all ([2007). For each realization T' c we 
can estimate T with a reasonable foreground template Tj 
which we suppose traces the morphology of the emission 
(see j gZQ) . 



Figure [9] shows histograms of 

T'a j = T(T>,T;)a j , 



(26) 



where aj is the standard deviation over unmasked pixels 
of the foreground template Tj (Tj = Ha, FDS, Haslam, 
and haze), for all 100 realizations of the CMB sky. We 
multiply by o~j to obtain the contamination in a "typi- 
cal" pixel in the ILC map. Each histogram has a mean 
consistent with zero, 

(r / ) « for each Tj (27) 

as expected since the CMB realization T' c is just as likely 
as — T' c . The implication is that we can only estimate the 
uncertainty in Tilc? we cannot explicitly correct for the 
bias. 

It is important to note that, although L cannot be 
corrected for the cross-correlation bias, the net effect is 
to decrease the variance relative to the true CMB, T c . 
Since the bias of the variance is not mean zero, the effect 
on the power spectru m of the CMB can be estimated (cf. 



on tne po wer spectrun 
iHinshaw et al.ll2QQ7[ ). 



4.2. HF type estimators 

The bias in a HF estimator is easy to understand and 
due entirely to the residual foreground emission after 
subtraction of thermal dust and free-free. Assuming per- 
fect subtraction of these two foregrounds, whatever other 
foregrounds are left in the HF estimator are subtracted 
off all lower frequency bands leading to a systematic bias 
towards softer spectra for the soft synchrotron and haze 
components. 

4.3. Forecasts for Planck 

The principle advantage of the Planck mission over 
WMAP is the large range of frequency coverage. In par- 
ticular, the multiple channels at very high frequency re- 
duce the foreground problem to (mostly) a single emis- 
sion mechanism, thermal dust. In order to exploit this 
feature, we suggest a new approach based on an ILC 
type CMB estimator. However, instead of minimizing 
the variance (which we have seen leads to large cross- 
correlation uncertainties), we choose the ILC coefficients 
(i to cancel as many power law foreground components 
as possible. This estimator can then be subtracted off of 
lower frequency bands with very minimal bias contami- 
nation of the inferred foreground spectra 

Planck's frequency coverage is 30, 44, 70, 100, 143, 217, 
353, 545, and 857 GHz. If we estimate that spinning dust 
emission is negligible above 100 GHz, then these bands 
contain mostly thermal dust (T oc with 1.6 < (3+ < 
2.3) and also small amounts of free-free, synchrotron, 
and haze (T ex v&- with -3.1 < /?_ < -2.1). 

In order to find an ILC that optimally cancels out this 
range of power law indices, we formulate the error func- 
tion, 

M/3) = (28) 

for a foreground spectrum, normalized to 94 GHz, T^ = 
T 94 (^/94 GHz)^, where Vi are the Planck bands above 
100 GHz and the sum is over bands. The physical inter- 
pretation of 894(13) is that it represents the fractional bias 
in the inferred 94 GHz amplitude of a foreground with 
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Fig. 8. — Integrated haze (from the bottom panel of Figure [4} in kJy/sr as a function of radial distance south of the Galactic center. 
The radial bins are 20 degrees wide and separated by 1 degree in longitude. The inner error bars are due to the formal error on the fit 
coefficients, the outer error bars are the 1-a standard deviation of the temperature fluctuations in a given radial bin, and the dotted lines 
represent the bias due to chance correlation between the CMB and the haze. 



power law index f3. We minimize J with respect 

to d over the range (3+ and /?_ defined above. Figure [T0l 
shows the resultant £94. With these Q we can form the 
Planck ILC Pp = CiPi, where Pi are the Planck data 
in thermodynamic mK at band i. We find that 



L P = -1.49Pi43 
+O.O2P545 - 



h 3.21P 217 - O.74P353 
3.21 x 10- 5 P 8; 



857 



(29) 



minimizes the area under the 5g 4 curve over the f3 range 
of interest. Using only the high frequency channels in 
this way avoids ambiguities due to an uncertain spinning 
dust spectrum. 

To estimate the contamination of Pp by synchrotron 
and haze emission we must make a couple of assump- 
tions. First, let us assume that these foregrounds follow a 
power law with with amplitudes of roughly our fit results 
at 23 GHz from Figure El and indices -3.1 < f3 s < -2.7 
and —2.7 < /3h < —2.4 respectively. Second, we make 
the implicit approximation that the morphologies of the 
foreground emission mechanisms do not change with fre- 
quency and that they are still well represented by the 
templates in §2.11 

Figure [10] shows our estimates for the error bars on 
the "mock true" soft synchrotron and haze spectra at 
WMAP frequencies given three CMB estimators: a 
WMAP only ILC (leading to cross correlation bias er- 



rors), our WMAP HF estimator CMB4, and the P P es- 
timator. As we have already seen, the WMAP only ILC 
clearly yields very large uncertainties, though they have 
the advantage of being mean zero (averaged over an en- 
semble of CMB realizations). CMB4 does significantly 
better, though now the errors are systematically biased 
towards softer spectra. With Pp, the bias errors are al- 
most completely eliminated. For a foreground with in- 
dices (3 = -3.1, -2.7, and -2.4, the bias errors at 23 (94) 
GHz are -0.11%, -0.06%, and +0.13% (-7.6%, -2.2%, 
and +3.3%) respectively. 

We emphasize that we have by no means attempted to 
formulate the "optimal" foreground removal algorithm 
for Planck. There have been numerous discussions on the 
topic describing and comparing different cleaning meth- 
ods. 13 Rather, we have demonstrated that, by exploiting 
the many bands and high frequency coverage, this sim- 



13 E.g., Brandt et al. (1994) undertook an earl y study of fore- 
groun d remova l with multiple frequency bands, Tegmark et al. 
( 2003) and Ide Oliveira-Costa fe Tegmarkj (|2006T ) describe an ILC 
type method in which the weighting depen ds on the multipo le £ 
of the spherical harmonic de compostion, Barreiro et al. ( 2004]) ex- 
pands on the MEM analysis, Eriksen et al. (2006) describe a com- 
ponent separation method, Hansen et al. (2006) tackle the problem 
with a wavelet analysis, and Dodelson (1997) attempts to "design" 
an experiment with optimal frequency coverage for removing fore- 
grounds. 
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Fig. 9. — Distribution of ILC chance cross correlation biases for 100 realizations of the CMB sky. a.) Hen-correlated emission, b.) 
dust-correlated emission, c.) Haslam-correlated emission, and d.) haze emission. 



pie Planck ILC will substantially reduce the systematic 
biases in the inferred foreground spectra compared to 
WMAP. 

5. DISCUSSION 

We have carried out a multi-linear regression fit to the 
WMAP data using foreground templates for free-free, 
soft synchrotron, and dust emission. We perform fits over 
both the (nearly) full-sky as well as smaller regions of 
interest. Our method simultaneously fits the amplitude 
and spectrum of each emission component and is immune 
to priors or initial "guesses" for the template amplitudes. 

Importantly, we find that the spectra of the foreground 
emission mechanisms cannot be determined to high ac- 
curacy with the WMAP data. The root cause is that any 
estimator of the CMB will necessarily be contaminated 
by the foregrounds to some degree, so that subtracting 
the CMB estimator from the WMAP data systemati- 
cally biases the inferred foreground spectra. Thus, the 
fit spectrum for each foreground is contaminated by some 
amount of CMB spectrum. The degree of systematic bias 
in the foreground spectra can be quite large and varies 
dramatically among our different CMB estimators. It is 
worth noting that the degree of contamination of a given 
CMB estimator does not significantly affect the variance 
of the estimator. 

We find that, upon removing the free-free, dust, and 
soft synchrotron emission, the 3-year WMAP data still 



contains the "haze" seen by Finkbeiner (2004). Since this 
haze emission is present in all residual maps regardless 
of which CMB estimator is used, we conclude that it 
is not an artifact of imperfect CMB subtraction. We 
have included a simple 1/r template for the haze to relax 
the stress on the other foreground components in our fit. 
This template imperfectly matches the morphology of the 
haze, but allows us to simultaneously fit an approximate 
amplitude and spectra of the haze emission as well. 

Despite the above issue of bias errors, we can make the 
following concrete statements about our fit foreground 
spectra. 

1. We find that the Ha-correlated emission does not 
follow a simple I v oc ^ _015 free-free powerlaw as 
expected from theory. Rather we find a bump in 
the Ha-correlated emission which we argue in our 
companion paper D+F (2007) can be explained by 
a mixture of free-free gas and spinning dust. We 
emphasize that this bump is not due to the contam- 
ination of the CMB estimator by the foregrounds. 

2. The dust-correlated emission has the now familiar 
fall off from 23 to 41 GHz consistent with a super- 
position of D+L spinning dust spectra and then a 
rise at 94 GHz from thermal dust emission. We 
cannot use our multi-region fits to map out varia- 
tions in the spinning dust spectrum from place to 




Fig. 10. — Left panel: fractional error in the inferred 94 GHz amplitude of a foreground with Toe ^ due to contamination of the Planck 
ILC Lp (see Equation I29|) . Center and right panels: biases due to foreground contamination of different CMB estimators for mock soft 
synchrotron (Tg oc i/Ps with —3.1 < (3s < —2.7) and the haze (Tjj oc v$ H with —2.7 < (3h < —2.4). The mock spectrum is shown with 
a solid line. The dotted lines are the mean zero ILC bias errors using only WMAP data alone, the dashed lines are the systematic biases 
due to CMB4, and the dot-dashed lines are the systematic bias using Lp. The systematic biases with Planck are more than an order of 
magnitude smaller than for a WMAP HF estimator. 

place across the sky since the uncertainty from the spectra of the foreground emission, a more reliable esti- 
systematic bias dominates. mate of the CMB sky is required. 



3. Our soft synchrotron spectra are very sensitively 
dependent on which CMB estimator is used. For 
example, with the WMAP team's published ILC, 
the soft synchrotron spectrum actually turns up 
at high frequencies, while other estimators show 
a steadily falling spectrum (with T v oc v$ where 
(3 ~ —3.0.) This ambiguity is entirely due to the 
contamination of the CMB estimator by soft syn- 
chrotron type emission. 

4. The haze spectrum is similarly very uncertain due 
to CMB estimator bias, however, for a given es- 
timator, the haze spectrum is always harder than 
the soft synchrotron. For reasons outlined in H'3A\ 
we suspect that the haze emission is due to hard 
synchrotron from a separate component of very en- 
ergetic electrons near the Galactic center. Possible 
sour ces include product s of dark matter annihila- 
tion ([Hooper et al.ll2QQ7l ) or a single energetic event 
in the last million yea rs, e.g. a gamma-ray burst 
(|Broderick et all I2QQ7D . 

Upcoming experiments like Planck will significantly 
reduce the problem of CMB estimator bias in the in- 
ferred foreground spectra due to the improved high fre- 
quency coverage. In particular, for frequencies v > 100 
GHz the Planck data will be virtually free of all fore- 
grounds except for thermal dust. Removing this one 
bright component with five bands is significantly less am- 
biguous than removing three or four bright components 
with five bands. Furthermore, what little free-free, soft 
synchrotron, and haze emission there is at these high 
frequencies can be roughly eliminated in a Planck CMB 
estimator through the appropriate choice of ILC coef- 
ficients. We estimate that the systematic uncertainties 
in soft synchrotron and haze spectra will be reduced by 
more than an order of magnitude compared to WMAP. 

Galactic foreground emission represents the most sig- 
nificant contaminant in determining cosmological param- 
eters. Ironically, the reverse statement is also true and in 
order to make progress towards determining the "true" 
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